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ABSTRACT 

Techniques of zero-temperature field theory that have found application in the anal- 
ysis of field theory at finite temperature are revisited. Specifically, several of the 
results that are discussed are relevant to the study of symmetry-changing phase 
transitions and high temperature QCD, which today are among the most actively 
investigated problems in finite temperature field theory. 

1. Overview 

Most of the techniques used in finite temperature field theory were first de- 
veloped for analysis of field theory at zero temperature. The purpose of these 
pedagogical lectures is to revisit (a part of) this field theoretical background. 

We start, in Sees. 2 and 3, by discussing techniques for studying the phase 
structure of a field theory using Green's functions generating functionals. These 
techniques are receiving renewed interest because of their role in numerous recent 
studies of temperature-induced symmetry-changing phase transitions, which have 
followed the suggestion that the observed baryon number asymmetry might have 
originated at the electroweak phase transition. 

Also relevant to the investigation of temperature-induced phase transitions is 
the formalism of non-equilibrium quantum field theories, which we discuss in Sec. 4. 

Finally, in Sec. 5 we review Chern-Simons theory, which has recently found 
application in finite temperature field theory because the Chern-Simons eikonal 
gives the generating functional of the hard thermal loops in QCD. 



2. Phase Structure and Generating Functionals 

2.1. Effective Action 

An important aspect of a physical system is its phase structure (what are the 
available phases, in what phase the system is at a given time, what is the nature 
of the phase transitions ....). Let us assume, at least at the beginning, that one 
of the fields in the theory ($) is a scalar and study the phase structure in terms 
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of the thermal expectation of that scalar field, defined by (we use units such that 

h = C = kBoltzman = l) 

tr e-P H § 

< <j> >= (1) 

tr e-? H ' 

where (3 is the inverse temperature 1/T. At zero temperature one is interested in 
the zero-temperature limit of < $ >, i.e. the vacuum expectation value < 0|$|0 >. If 
one could solve the theory completely these expectation values could be evaluated 
explicitly from their definitions. Unfortunately there are only a few quantum field 
theories that can be solved completely (and they are not very interesting from the 
physical point of view ), so that the answer to the problem at hand can only be 
approached somewhat indirectly by setting up additional formalism, with the hope 
that a resonable approximation scheme emerges. 

We use an approach based on generating functionals: the effective action, the ef- 
fective energy, and the effective potential. This was first proposed by Jona-Lasinio 1 , 
and later developed in Refs.[2-5]. (For a review of this subject see Ref.[6].) For 
definiteness, we work at zero-temperature; however, the finite temperature general- 
izations can be easily obtained following the discussion in Sec. 3. 

The first step is the introduction of the "partition function" Z{J) 



i f J(x)$(x)\0 > , (2) 

J X 



Z(J) =<0|Texpi / J(x)<f>(x)\0 > 

where T is the time-ordering operator. 

Z(J) is a derived quantity, a quantity external to the theory, but its importance 
is due to the fact that Z(J) is the generating functional for Green's functions: 



G n (x u ,x n ) =< 0\T<f>( Xl ) $(x n )\0 >= 



5 n 

Z(J) 



i n 5J{x 1 )....J(x n ) 



(3) 

,7=0 



In as much as all the physical information on a quantum field theory is encoded in 
the Green's functions, one can study the theory just using the generating functional 
Z(J). 

The generating functional for connected Green's functions W{J) is given by 

W(J) = -ilnZ(J) . (4) 

The same physical information encoded in the Green's functions is organized more 
economically in the connected Green's functions; in fact, once the connected ones 
are known the complete Green's functions can be reconstructed. 

The effective action r(0) is defined as the functional Legendre transform of W(J): 

r(0) = w(j) - f 4>{x)j(x) , (5) 

J X 



where 



(6) 



and in (5) J is expressed in terms of <j> using (6). [T(4>) is sometimes referred to as 
the "quantum action" ; in the limit K — > it reproduces the classical action.] 
Evidently 

-'<*>• « 7 » 

and, as a consequence, physical solutions, which require vanishing J, correspond to 
the stationary points of r. 

From the definitions that we have introduced it is clear that 

<0|$(x)|0 >= [<£(*)] j=o = Mx), (8) 

and therefore using Eq.(7) one finds that the vacuum expectation value of $ is given 
by the value of <j)(x) which stationarizes T ((/>): 



5<j>(x) 



= 0->Mx) =<0|$(x)|0> . (9) 

4>(x)=4> (x) 



We have replaced the problem of evaluating the vacuum expectation value of a 
quantum field by the problem of finding the stationary points of a classical func- 
tional. Such a problem has a long tradition in physics and we can use the numerous 
approximation techniques that have been developed over the years. 

Finally we note that T(4>) is the generating functional of the one particle irre- 
ducible* (1PI) Green's functions: 



(5 n r((/>) 



S<l>(xi) S(j>(x n ) 



= Gl l PI (x 1 , ,x n ) . (10) 

4>(x)=4> (x) 



These contain all the physical information encoded in the Green's functions in an 
even more economical format than the connected Green's functions. 

2.2. Effective Potential and Effective Energy 

If one is interested in all the physical information in a quantum field theory, it 
is necessary to study the complete effective action, as a functional of the arbitrary 
background field 4>{x). However, when one is only interested in the vacuum expec- 
tation value 4> (x) =< 0|$(x)|0 >, one can use any a priori information available on 
(f>o(x) to simplify the analysis. 

If it is expected that translation invariance is not broken, as is usually the case, 
one can assume that <t>o(x) be independent of the space-time point: <j>o{x) = constant = 
fa. With this hypothesis, in order to determine fa it is sufficient to study the effective 
action as a function of a constant background fa However, when <f> is constant an 
infinite volume factor arises from the space-time integrations, and therefore, rather 
than working with the effective action, we introduce the effective potential V(fa, 
defined by 

V(fa [ = -[T(^)}^ constant . (11) 



A graph is said to be "one particle irreducible" if it does not become disconnected upon opening any 
one line; otherwise it is "one particle reducible" . 



The effective potential is the generating functional of the 1PI Green's functions 
at zero energy and momentum*, and in the h — > limit the effective potential 
reproduces the classical potential. Clearly from Eqs.(9) and (11) it follows that for 
translationally invariant theories the vacuum expectation value can be determined 
using the effective potential: 

= O^0 o =<O|$(a;)|O> . (12) 

= 00 

An intermediate possibility between evaluating the effective action as a func- 
tional of a space-time dependent background <f>(x) and the much simpler case of 
a constant background, is the evaluation of the effective action for static (time- 
independent, position-dependent) background </>(x). In this circumstance, the time 
integration in r(0) leads to an overall (infinite) time-interval factor; it is therefore 
convenient to introduce the effective energy E(<f>), defined by 

E(<j>) J dx a = -[r(0)] 0=stQtJC . (13) 

The effective energy is the generating functional for 1PI Green's functions at zero 
energy, in the h — > limit it reproduces the classical energy, and in time translation- 
invariant theories it can be used to determine the vacuum expectation value using: 

= O^0 o (x) =<0|$(a;)|0> . (14) 

0(x)=0 O (x) 



5V{4>) 
5(f> 



5E{4>) 



2.3. Loop Expansion 

In the search for the stationary points of T{4>), the choice of a strategy for eval- 
uating T(4>) is very important. One commonly used technique is the approximation 
of T{4>) based on the following known 3 * 4 expansion** 

r(0) = I(0)-ltrlnD + r 2 (0) . (15) 
I{4>) is the classical action. The functional operator D^ 1 ((/)■, x, y) is defined by: 

-TOT < 16 » 



Note that in general the Green's functions at zero energy and momentum are not physical, though they 
might have physical relevance in particular contexts. The poles and the residues at the poles of the Green's 
functions are physical. 

** The expansion is formally in powers of fi but this doesn't explicitly appear in Eq.(15) because in our 
units h = 1. 



The term |[tr In £>], where the trace and the logarithm are taken in the functional 
sense, is the one loop contribution [it is 0(h)]. T 2 (<f>) is computed as follows. In the 
classical action /(<&), shift the quantum field by a c-number field 4>{x). /($ + </>) 
contains terms cubic and higher in $ that define an interaction part I int (^;^) whose 
vertices depend on 4>. T 2 {4>) is the sum of all the two-or-more-loop 1PI vacuum 
graphs in the theory with vertices given by I in t(4>',^) and propagators taken to be 
D(cp;x,y). 

For example, for the A$ 4 scalar theory 



/(*) 



4! 



D-\fr x,y) = (-□ + a - ^ 2 (x))5 4 (x - y) , 



hnt(4> 



;*)= / 

J X 



_^^$3 _ A 4 

6 V 4! 



(17) 

(18) 
(19) 



The diagrams which contribute to T 2 (4>) for the A$ 4 scalar theory are shown in Fig.l 
up to three loops. The lines represent D((f>;x,y) and, as indicated by Eq.(19), there 
are three-point and four-point vertices. 
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e 

Fig. 1 

A loop expansion for the effective potential can obviously be obtained by eval- 
uating Eq.(15) for constant and extracting a (infinite) volume factor: 

V(4>) = V cl {4>)+ l -J In D{4>; k) + V 2 (<P) , (20) 

where V c i(4>) is the classical potential, £>(</>; k) is the Fourier-transform of D(<fr,x,y) at 
constant <j) 

D(frx,y)^ f e-^-y)D(frk) , (21) 

Jk 

and V 2 (4>) is given by 

V 2 (<f>) / = -[T 2 (<f>)] ^constant ■ (22) 

J X 

Analogously, the loop expansion for the effective energy is obtained by evaluating 
Eq.(15) for static 4> and extracting a (infinite) time-interval factor. 





2.4- Selective Summations and Effective Action for Composites 

It is very difficult to carry out the loop expansion beyond two loops. It is 
however possible to perform selective summations of higher loop graphs. One way 
to systematize such summations is the large TV method, in which one considers an 
iV-component field and uses the fact that in the large N limit some multi-loop graphs 
give greater contributions than others. For example in the A<I> 4 TV-component scalar 
theory the leading multi-loop contributions come from bubble-graphs (see Fig. 2). 



Fig. 2 

Another way to perform systematic selective summations makes use of the effec- 
tive action for composites 7 . This generalization of the effective action was primarily 
introduced to seek the phase structure at the level of the two-point function*, in 
addition to, or instead of, the one-point function. Moreover, it can also be used to 
rexpress the ordinary effective action in a way that is more suitable for selective 
summations. 

One arrives at the definition of the effective action for composites following a 
path very analogous to the one that leads to the ordinary effective action. The 
starting point is a new partition function that depends on a bilocal source K(x,y), 
in addition to the local source J(x): 



Z(J,K) = e lW ^ =< 0\Texpi 



f J(x)$(x) + l f $(x)K(x,y)$(y) 

Jx ^ J x,y 



|0 > , (23) 



where we have also introduced W(J, K). A possible expectation value of the quantum 
field $(ar) is <j>(x); this is again defined as the functional derivative of W with respect 

but now one also introduces G(x,y), defined by 



For example Fermi fields do not have a vacuum expectation value, but the fermion propagator can show 
spontaneous chiral symmetry breaking through mass generation. 



which is a possible expectation value of the connected part of the two point function. 

The effective action for composites T(<j>, G) is defined as the double Legendre 
transform of W(J, K) 



T(<f>,G) = W(J,K) 



4>{x)J(x) 



(p(x)K(x,y)(j>(y) 



x,y 



G(x,y)K(y,x) 



(26) 



x.y 



where J and K are to be expressed in terms of 4> and G using (24) and (25). 

One can show 7 that T(<j>, G) is the generating functional in <f> for two-particle 
irreducible* (2PI) Green's functions, expressed in terms of the full propagator G. 

From Eq.(26) it is evident that 



54>(x) 



= -J{x)- f K(x,y)<Kv) , 



SG(x,y) 2 K{X ' V >- 



(27) 



(28) 



Since physical processes correspond to vanishing sources J and K, Eqns.(27) and 
(28) imply that T{4>, G) is physical at its stationary points. 

Eq.(28) also indicates the relation between T{<j>,G) and the ordinary effective 
action; in fact, as it is clear from its definition, T(<cj>, G) must reproduce T{4>) when 
K = 0, and using Eq.(28) one finds that 



r(0) = r(^,G o ) 
when G satisfies the stationary requirement 



. 



.SG(x, y) j g=Gq 

A loop expansion is also available 7 for the effective action for composites 

T(<f>, G) = Id(<f>) - % -tr[\nG - D~ l G] + T 2 (cf>, G) , 



(29) 



(30) 



(31) 



where r 2 (0, G) is given by the sum of the two-or- more- loop 2PI vacuum graphs in 
the theory with vertices given by I int {4>;§) and propagators set equal to G(x,y). 

For G = G Eq.(31) gives a new expansion for the ordinary effective action. 
Moreover, using Eq.(31) the gap equation (30) which determines G can also be 
expressed in terms of a loop expansion as follows 



G \x,y) =D- 1 (x,y)-2i 



SG(x,y) 



(32) 



G=G 



A graph is said to be "two particle irreducible" if it does not become disconnected upon opening any 
two lines; otherwise it is "two particle reducible" . 



Eq.(32) gives the relation between G and the tree-level propagator D, and can there- 
fore be used to express the diagrams of the expansion (31) in terms of diagrams of 
the ordinary expansion (15). For example, for the A<1> 4 scalar theory it is easy to show 
that the summation of all the contributions of the bubble graphs to the ordinary 
effective action r(</>) corresponds to the solution of the Eqs.(29)-(30) when r 2 (<f>, G) is 
approximated by its only graph of 0(A), i.e. the double-bubble graph in Fig.l. This 
bubble-resummation has been used in a recent study of the temperature-induced 
phase transition of the A$ 4 scalar theory 8 . In general the solution of Eqs.(29)-(30) 
with the approximation of r 2 (<j>, G) by a finite sum of its graphs corresponds to the 
selective summation of an infinite class of "ordinary graphs" [graphs of the ordinary 
expansion (15)]. 

When interested in translation-invariant (time translation-invariant) solutions 
one can perform selective summations of multi-loop graphs which contribute to the 
effective potential (energy). These summations may be performed systematically 
using the effective potential for composites V(<f>, G) and the effective energy for com- 
posites E((p,G), which, in analogy with Eqs.(ll)-(13), are defined in terms of the 
effective action for composites by 



2.5. Variational Interpretation and Variational Approximations 

We have seen that, in the formalism of these generating functionals, physical so- 
lutions satisfy variational equations which stationarize the functionals. In this sub- 
section we shall discuss the fact that the functionals themselves have a variational 
interpretation. The knowledge of these variational interpretations is important, 
especially because they can suggest approximations. 

Let us start with the effective energy E (</>). It can be shown that E(4>) is the 
minimal expectation value of the Hamiltonian in normalized states which are con- 
strained so that the expectation value of the quantum field $(ar) in those states is 
equal to </>(x), i.e. 



This result can be immediately related to the general quantum mechanical prin- 
ciple that in searching for the solutions of a static problem one should minimize 
the expectation value of the Hamiltonian. Using the effective energy this princi- 
ple is implemented in two steps: one first minimizes the expectation value of the 
Hamiltonian with the constraint < *|$|* >= <f>, and then the resulting ^-dependent 
functional (the effective energy) is minimized with respect to <f>. 

For the effective potential V{4>) a similar variational interpretation and similar 




(33) 




(34) 



E{4>) =min < > , 

< W|W >= 1 , < *|$(a:)|* >= 0(x) . 



(35) 



considerations apply; one finds that 



V{<j>) 
< >= l 



mm 



< *|$(a;)|* > 



(36) 



Also the effective potential and the effective energy for composites have a varia- 
tional interpretation. For example to obtain E(4>,G) one minimizes the expectation 
value of the Hamiltonian on normalized states which satisfy two constraints: 



E(<j>, G) =min < ^\H\^ > , 
< >= 1 , < >= <f>(x) , 

< >| Xo=yo = 0(x)0(y) +G(x,y) . 



(37) 



The variational definitions for the effective potential and the effective energy are 
related to the general quantum mechanical variational principle for static problems. 
Correspondingly, the variational definition for the effective action T(4>) can be derived 
from the Dirac variational principle 9 for time-dependent problems, which instructs 
minimizing the following quantity 



J dt <y,t\id t - H\V,t>= J dtdx 



**(t,x)(idt-H)y(t,x) 



with the normalization constraint 



<tf,t|#,i>= / dx#*(i,x)tf(t,x) = 1 . 



r(<£) can be variationally defined as follows 



J dt < $-,t\id t - H\V+,t> 



(38) 



(39) 



(40) 



where the states |*+,t > and > are constrained so that 

<*_,t|$(x)|*+,t> = 0(t,x) 
< #_,t|* + ,i > = 1 
lim |tf±,i > = |0 > . 

t — s-±00 



(41) 



The variational definition of the effective action for composites T((j>, G) is analo- 
gous, but with the addition of a constraint for < vE/_ 7 <|$( x )$(y)|* + ,t >. 

The variational interpretation of our generating functionals naturally leads to 
some variational approximations. These approximations are best formulated in 
the Schrodinger picture for quantum field theory 1011 , in which states are described 
by wave functionals of a c- number field 4>(r), the inner product is defined by 
functional integrals 

< ^1*2 >~ f r>0tfi(<£)tf 2 (0) , (42) 



the action of the field on the state is by multiplication 



>~ 4>^{<f>) , (43) 

and the action of the canonical momentum on the state is by functional differenti- 
ation 

n|* >~ \j^{<t>) • (44) 

In this formalism one finds the nice result that if ^ G {4>) is a Gaussian functional of 
the form 



#g(</>) = exp 



(45) 



-\Jdx dy (0(x) - 0o(x)) G-^x.y) (0(y) - 0o(y)) 
then the expectation of the Hamiltonian in *g(0) automatically gives 

< H > G = J D4®* G (<I>)H* G (<I>) = E bubble (cf> , G) , (46) 

where E bubble {<t> ,G) is the approximation of the effective energy for composites which 
includes only the 0(A) contribution to E 2 (<j> ,G). Therefore by stationarizing < H > G 
with respect to variations of G one obtains the "bubble-resummed" effective energy. 

2.6. Comments and Cautions 

Some caution is necessary in using the generating functionals that we have dis- 
cussed. Let us start by analyzing some features of the effective potential* in the 
simple A$ 4 scalar theory. Including only the classical and the one loop term one has 
(see Sec.2.3) 

V{<t>) ~ V cl {4>) + \J k ln[fc 2 - a + ^ 2 ] , (47) 

where we have set up a Euclidean evaluation of the one loop contribution. (N.B. We 
are not doing Euclidean field theory. We simply evaluate a Minkowski space-time 
integral using a continuation to Euclidean space.) 

The first important observation is that the integral in Eq.(47) is clearly diver- 
gent, and therefore a regularization/renormalization procedure is needed. It turns 
out that this renormalization can be performed in a rather standard fashion 2-5 . 

Eq.(47) also indicates that the effective potential is not necessarily real and con- 
vex (if a > the effective potential in Eq.(47) is not convex and has an imaginary part 
for small 4>) , and therefore it cannot be interpreted as the free energy of the system. 
It has been suggested 12 that the imaginary part of the effective potential might carry 
information on the decay properties of the system. The possible non-convexity of 
the effective potential is important in the Legendre transform procedure 13 ; in fact, 
it may not be possible to find a unique solution for J(4>) in Eq.(6), because, e.g., 
there are multiple roots. In these cases, if one is interested in the free energy, a 



Clearly the issues discussed in this section also concern the effective action and the effective energy, but 
for definiteness we base the presentation on the case of the effective potential. 



Maxwell construction can be used, but then the resulting effective potential will not 
be the generating functional for the Green's functions. 

Concerning the effective potentials of gauge theories there is obviously an ad- 
ditional issue to be addressed: gauge invariance. In a gauge theory the effective 
potential, besides depending on a background scalar field (and its complex conju- 
gate), can also depend on a background gauge field: V(A li , 4>,4>*). There are two ways 
in which the effective potential may fail to be gauge invariant. First it may not be 
invariant against gauge transformations of A^, <fi, 4>*. Second, even if the effective 
potential is invariant against gauge transformations of A M , 0*, it may still depend 
on the gauge choice made in quantization; for example in covariant gauges it may 
depend on the coefficient £ of the gauge fixing term ^(d^A^) 2 . In general the effective 
potential does depend on £. 

The gauge dependence of the effective potential can be understood in terms of 
the fact that off-shell Green's functions are not physical and are not gauge invariant. 
Although the full effective potential need not be a physical quantity in a gauge 
theory, gauge invariant physical information can be obtained from it, for example 
at the stationary points. However, in many cosmological models the detailed profile 
of the effective potential, even away from the stationary points, plays an important 
role (slow rolling down out of equilibrium). Therefore for these models the physical 
interpretation of the full effective potential is an important issue. 



3. Finite Temperature Relativistic Quantum Field Theory 

The observable properties of a relativistic quantum field theory at nonzero tem- 
perature* can be extracted from the thermal Green's functions, defined by 

Gp(x u ,x n ) =< T<f>( Xl ) $(a; n ) >= / — - . (48) 

tr e 

Let us consider the two-point thermal Green's function G@(xi,x 2 ). The Dyson- 
Schwinger equations of motion for Gp(xi,X2) follow from the Heisenberg equations 
of motion like in the ordinary zero temperature case 

n 1 G (x 1 ,x 2 ) = -i5( Xl - x 2 )+ < rn$(xi)$(a; 2 ) > . (49) 

The structure of Eq. (49) is exactly the same as in the zero temperature case. The 
difference between the thermal Green's function and the zero temperature Green's 
function is merely in the boundary conditions, which determine the unique, phys- 
ically appropriate solution to the differential equation. One can therefore follow 
the same method of solution to Eq.(49) used at zero temperature, except for what 
concerns the boundary conditions. 

For definiteness and simplicity, we illustrate the methods of solution of Eq.(49) 
in the case of bosonic fields, for which the boundary conditions are 

Gp(x-L,x 2 )\ x o =0 = Gp(x 1 ,x 2 )\ x o=_ i0 . (50) 



* For a review, see Refs.[14, 15]. 



These follow from the cyclicity of trace, which defines Gp, see Eq.(48). [Note that 
Eq.(49) is in real time, but the solution is required to be periodic in imaginary time.] 
We also assume that, as it happens in most cases of physical interest, translation 
invariance is not broken, and therefore 



Gp(x 1 ,x 2 ) = Gp{xi - £2,0) = Gp{xi - x 2 ) ■ (51) 

Because of the translation invariance it is convenient to go to conjugate space. 

In the "real time formalism" one introduces the Fourier integral transform Gp(p) 
of Gp(x! - x 2 ): 

Gp( Xl -x 2 )= f e-W^-^Gplp) . (52) 

It then follows simply from the fact that the field $(x) is an operator which verifies 
causal commutation relations, that Gp{p) has the following integral representation 

Gp{p) =i r ^-^V + /(^(Po,P) 

(53) 



=>F 



d P0 (I \ 

(2^)^(Po.P) 



l + fip'o) f(Po) 



Po-p' + ie Po-Pq- it 



where / is the bosonic distribution function 

and pp is a density function to be determined using Eq.(49). 

In the lowest order in perturbation theory the density function is known, 

p f / ee = 27Te(p )S(p 2 - m 2 ) , (55) 

and correspondingly the thermal 2-point Green's function is given by 

G f'**< rt ^7^ T - + ^* 2 - m J ). (50) 

[Here e(po) = 0(po) - Q(—po), and 9 is the unit step function, vanishing for negative 
argument.] Note that G^ ree satisfies the free- field Schwinger-Dyson equation (p 2 - 
m 2 )G free = i, and the temperature dependence is confined in the term, proportional 
to S(p 2 - m 2 ), in which boundary conditions are encoded. 

As an alternative to the real time formalism one can work in the "imaginary 
time formalism" , in which the Green's functions are considered in imaginary time. 
As a consequence of the periodicity in imaginary time described in Eq.(50), in 
this formalism the momentum space two-point thermal Green's functions are the 
Gp(u n ,p) defined by Fourier series and integral transforms 

[G (x)] Xo =- iX4 = ^E e ^ 4 / (^3 eipX ^K'P) > (57) 



where 

u„ = i— . (58) 

In Eq.(57) there is a Fourier series for the time-like component and an ordinary 
Fourier integral transform for the space-like components. 

The same analysis which leads to the spectral representation of the real time 
Green's functions also leads to the following spectral representation for the imagi- 
nary time Green's functions: 



Gp(u n ,p) = i I 
J—t 



(2tt) u n - p' 



(59) 



which involves the same spectral function pp present in Eq.(53), and the same 
denominator structure, but now p is substituted by the imaginary discrete quantity 
u n and therefore, because the denominator never vanishes, there is no need for ie- 
prescriptions. 

The imaginary time Green's function which corresponds to the spectral function 

P fr is 

ffiee K> p ) = __J__ , (60 ) 

which looks like a free zero temperature Green's function, except that p has been 
replaced by u n - 

Perturbation theory has its most physical realization in the real time formalism, 
however some technical problems are encountered. The structure of perturbation 
theory is the same as at zero temperature because the structure of the equations 
is the same, but it is difficult to implement the boundary conditions. Moreover, 
the Green's functions involve a <5-function and therefore in high order calculations 
one encounters ambiguous products of ^-functions at the same point. (Only the 
one-loop results can be obtained rather easily in the real time formalism because no 
products of ^-functions at the same point can occur.) Obviously these difficulties 
can be overcome by doing things carefully, but it is very hard to find algorithms 
for perturbative calculations in the real time formalism. A systematic algorithmic 
prescription is given by the time path approach 14 , but involves mathematical elabora- 
tions that obscure physical content: a specific closed time contour in the complex 
plane is chosen, and one uses both time- and anti-time-ordered Green's functions 
with an effective doubling of the degrees of freedom. 

Perturbation theory is not pathologic in the imaginary time formalism. In fact, 
one can introduce Feynman rules in complete analogy to zero temperature, except 
for the substitution in the propagators of the continuous variable p a by the discrete 
u n , and the substitution of the integration over a loop four- momentum by an in- 
tegration over the three spatial components of the momentum and a sum over the 
discrete values of the time-like component: 



dp 

(2^)3 



Unfortunately the simplicity of perturbation theory in the imaginary time for- 
malism is somewhat offset by the fact that, when one is interested in the answers to 
dynamical questions, it is necessary to continue results back to real time. However, 
the effective potential, which is the generating functional for Green's functions at 
zero momentum, can be computed in the imaginary time formalism because at zero 
momentum (static physics) there is no difference between real time and imaginary 
time. For example for the A$ 4 scalar theory both in imaginary time and in real time 
one finds that the one-loop effective potential is given by 

= W) + \ £ / (^3 + k 2 - a + ^) . (62) 



4. Non-Equilibrium Quantum Field Theory 

4-1. Isoentropic Energy-non- conserving Time Evolution 
A physical system is described by its density matrix p, 

e -0H 

P^^r-m^ (63) 
and average values of observables O are determined by the density matrix: 

< O >= trpO . (64) 

Finite temperature field theory describes equilibrium physics, and therefore in- 
volves a time-independent density matrix. In general, however, the density matrix 
is time-dependent, and the task of non-equilibrium quantum field theory is to study 
the time evolution of p. 

Non-equilibrium physics is a vast subject, and there is no canonical approach to 
its investigation. Usually the approach is suggested by the specifics of the physical 
system that one wants to describe. We discuss an approach which is useful in early 
universe cosmology, and is set up in the framework of the field theoretic Schrodinger 
picture, which is particularly suitable to time-dependent problems that require an 
initial condition for specific solution. 

The functional density matrix is given by a superposition of wave functionals 

p{<h,<h) = 5>n*n(<M*;(0 2 ) , (65) 

n 

were {*„} is a complete set of wave functionals, and p„ is the probability (J2 n p n = l) 
that the system is in the state \I>„. In general both and p n are time-dependent. 

In equilibrium the dynamics is time-translation invariant and energy is con- 
served. The complete set of wave functionals {*„} can be chosen to be the set of 
the (time-dependent) energy eigenstates, and the p n are time-independent and are 
given by the canonical Boltzmann distribution: 

e -l3E n 

P" = v e-0E n ' (66) 

Z—in 



where E n is the energy eigenvalue of the state *„. The time evolution of the density 
matrix is trivial: it remains constant in time because both the p n 's and are 
constant (N.B. the time dependence of the *„'s is just a phase). 

For non-equilibrium physics the time evolution of the density matrix is instead 
nontrivial. In fact, the p„'s need not be Boltzmann factors and can change in 
time, and it might not be possible to choose the *„'s as energy eigenstates. We 
assume that the time dependence of the *„'s be determined by the time-dependent 
Schrodinger equation. As a consequence, the density matrix p satisfies the following 
differential equation 

n n n 

In order for Eq.(67) to describe a well defined initial value problem for the time 
evolution of p, it is necessary to give the form of H and a model for dp n /dt. We in- 
vestigate Eq.(67) in the case of time independent p„'s (which corresponds to entropy- 
conserving time evolution) and time dependent Hamiltonian* (which corresponds 
to energy non-conserving time evolution) 16 . With these hypotheses Eq.(67) takes 
the form of the quantum Liouville-von Neumann equation with time dependent 
Hamiltonian 

f t =i[p,H]. (68) 

In particular, we assume that the "interesting" time dependence of H occurs 
in an interval U < t < t s whereas H = Hi = constant when t < u and H = H f = 
constant when t > tf. Interesting Hamiltonians, which may be investigated, include 
Hamiltonians with time-dependent mass squared (in particular the mass squared 
changing sign is useful in the study of symmetry-changing phase transitions), and 
Hamiltonians corresponding to a quantum field theory in a background Robertson- 
Walker metric with time-dependent scale factor. 

We impose as initial condition that the density matrix p when t < U be given by 
the Boltzmann distribution for Hi with temperature 7*. The typical result, which 
one seeks using the formalism that we developed, is the density matrix at times 
t > t f: i.e. one solves Eq.(68) and examines p at late times. If p still changes with 
time when t > tf, one concludes that the system remains out of equilibrium. If p 
is time-independent for t > tf, there are two possibilities: either p is given by a 
Boltzmann distribution for Hf with a temperature Tf, in which case the system 
regains thermal equilibrium at temperature Tf, or p is not given by a Boltzmann 
distribution, in which case the system reaches non-thermal equilibrium. 

4-2. Methods for Solving the Liouville-von Neumann equation 

We are interested in obtaining solutions to the Liouville-von Neumann equation 
(68) with a canonical density matrix as initial condition. However, the exact solution 

si. 

The physical idea behind these assumptions is that one drops the time dependence of the p n 's, in order 
to make the equations treatable, but tries to account for (part of) this time dependence by introducing an 
"effective" time-dependent Hamiltonian. 



of this problem is only feasible 16 for problems that are described by a quadratic 
Hamiltonian. For more general Hamiltonians one can obtain approximate solutions 
using the observation that the Liouville-von Neumann equation can be derived 
from a variational principle introduced by Balian and Veneroni 17 . An approximate 
application of this variational principle with a restricted variational Ansatz leads to 
approximate equations for the density matrix that may be integrated. 

The Liouville-von Neumann equation can be derived by varying the actionlike 
quantity 



dt tr 



[MM)] 



t=ti , 



(69) 



where A, p, and H are time-dependent kernels and the trace is over these kernels. A 
is a Lagrangian multiplier kernel. By varying I with respect to A and p one obtains 
respectively 



5aI 



6 P I 



-r 



dt tr 



dp 
~di 

dA 
~dt 



+ i[H,p] 
+ i[H,A] 



SA - [tr{ P 5A)] t=tf 



Sp - [tr(Adp)] 



t=u 



(70) 



(71) 



As in all time dependent variational principles we impose boundary conditions; we 
require that 

[SA] t=tf = , (72) 



so we may set 

[A] t=t/ = identity 

Also we require, according to our program, 



(73) 



[p]t=ti = Boltzmann distribution 



(74) 



and therefore 



[Sp]t=u = 



(75) 



Demanding that I be stationary against both variations of A and p, with the above 
boundary conditions, gives the Liouville-von Neumann equation for p, and also for 
A. The boundary condition Eq.(73) selects the static solution A = identity, for all 
time, so that A disappears from the discussion and we are left with a variational 
formulation of the Liouville-von Neumann equation for p. 

If, rather than performing arbitrary A and p variations, we evaluate / with spe- 
cific parameter dependent expressions for A and p and then vary these parameters, 
we obtain an approximate solution of the problem. In particular, a Gaussian Ansatz 
can be very useful. In this case p can be parametrized, also exploiting the Hermitic- 
ity of p [p(4>i,4> 2 ) = p*{4>2-,4>i% as 



p(^i,</>2) = e ^exp 



G~ 



2m 



G 



/2 



, (76) 



where G and n are real and symmetric kernels, £ is Hermitian, and 7 is to be chosen 
so that p is properly normalized. £ is called the degree of mixing; in fact, it is a measure 
of the amount by which p differs from a pure state. For £ = 0, p = *((/>i)**(</> 2 ), with 




- 2ffl U 



(77) 



The meaning of G and n can be understood by looking at the following bilinear 
expectation values 



This approximation scheme has been used in an investigation of early universe 
evolution 16 . The results support the validity of the heuristic method of analysis 
of early universe phase transitions (a dynamical problem) in which one evaluates 
only the (static) effective potential and assumes that the evolution be adequately 
described by substituting the effective potential in place of the classical potential 
in classical evolution equations 18 . 

5. C hern- Simons Theory 

5.1. Physical Relevance of 3- dimensional Quantum Field Theories 

In this section we discuss a particular 3-dimensional quantum field theory. There 
are several physical circumstances in which a 3-dimensional quantum field theory 
can be relevant. 

A first example is given by Minkowski 4-dimensional field theory in the Schrodin- 
ger Hamiltonian formalism, wherein one works at fixed time and therefore the c- 
number arguments of the wave functionals depend only on the three spatial coor- 
dinates, they are functionals of 3-dimensional fields. 

Another example is given by situations in which the interesting physics is (at 
least approximately) confined to a plane. This scenario occurs in some condensed 
matter phenomena, like the quantum Hall effect and high T c superconductivity, and 
for motion in presence of cosmic strings, which is described by planar gravity. 

Yet another circumstance in which 3-dimensional field theories are relevant is 
directly related to the topic of this workshop. In fact, one expects that in the 
high temperature limit thermal field theory in four dimensions can be described 
by a 3-dimensional field theory. A naive formal argument, which is often used to 
support this expectation, is based on the observation that, because of the periodic 
boundary conditions imposed on the fields, finite temperature field theory in the 
imaginary time formalism "lives" in R 3 x S 1 , and the radius of the S 1 is given by 
1/T. In the limit T — > 00 the S 1 collapses leaving a 3-dimensional space R 3 . This 
apparent dimensional reduction in the high temperature limit can be also seen in 



< n(r)n(r') > 



< $(r)Ii(r') > 



< $(r)*(r') > 



[ G i/2(i_£)-i G V2] (r y) 

i[G- 1/2 (l+£)G- 1/2 ](r, r') + 4[nG 1/2 (l-0 _1 G 1/2 n](r, r') 
^(r-r') + 2[G 1 / 2 (l-0- 1 G 1/2 n](r,r') . 



Feynman diagrams; for example in the A$ 4 scalar theory the tadpole contribution 
to the self-energy is 



tT ^I (2tJ)3 A -4n 2 7r 2 T 2 - p 2 - m 2 ' (?9) 

and its high temperature limit, in which only the n=0 term of the sum survives, is 
given by 

* "^3. (80) 



(2tt)3 p 2 , ,,. 

The passage from Eq.(79) to Eq.(80) can be interpreted again as a dimensional 
reduction (and a dimensionful scaling of the coupling A — > AT). 

This possible thermal field theory application was one of the initial motivations 
for studying the specific 3-dimensional field theory, involving the Chern-Simons 
term, which we review in the following. Recent results 19 in the context of the hard 
loops high temperature approximation have found that a quantity which naturally 
appears in the context of Chern-Simons theories is indeed directly related to the 
generating functional for hard thermal loops in QCD. 

5.2. The Chern-Simons Term 

The Chern-Simons term* Q(A) is a structure which can be constructed out of 
gauge fields in three dimensions: 

n(A) ee -J-^e^tridiAjAk + ^AA.Ak) . (81) 

(If the gauge fields are Abelian the last term on the right is absent.) Q(A) is a 
topological quantity, it is independent of the metric. 

Even though it is a 3-dimensional object, the Chern-Simons term appears in 
the study of some 4-dimensional problems. It arose first in the Schrodinger-picture 
description of the vacuum angle of quantized 4-dimensional non- Abelian gauge the- 
ories. This comes about in the following manner. Consider the quantity 

W(A) = J dx Q(A) , (82) 

which satisfies Gauss law, i.e. 

°^X w{A)=Di (^ eijkF ? k ) =o - m 

However W(A) is not gauge invariant: under a finite gauge transformation, corre- 
sponding to the gauge group element g, W(A) changes by the winding number n g of 
the gauge transformation: 

W(A") - W(A) = J rfx [Q(A") - Q(A)} - J dx JLe^c^ g' 1 A k ) + u{g)\ = n g , (84) 



For a review, see Ref.[20]. 



where 

"GO = ^^{g^dig g~ l d 3 g g- l d k g) , (85) 

n g = J dx u(g), and we assume that the gauge field vanishes sufficiently rapidly as 
|x| — > oo so that the surface term, which arises from / dx e i3k di(djgg~ 1 A k ) : also vanishes. 
It is also true that u(g) can be written as a total derivative, i.e. there exist some 
uj i (g(q)), where q(x) is a parametrization for g(x), such that 

io(g) = d t io\g(q)) . (86) 

However, because of the specific behavior of the parametrization q of the group 
elements as |x| — > oo, the integral of over the whole volume is not zero, in spite 
of the fact that, owing to Eq.(86), it can be presented as an integral over the surface 
at infinity. Rather the integral takes an integer value n g , which characterizes the 
homotopic equivalence class to which g belongs. 

The relevance of all this to 4-dimensional QCD can be seen in the Schrodinger 
representation, where states are 3-dimensional functionals that must satisfy Gauss 
law. Therefore they can take the form 

= e tew{A ^{A) , (87) 

where ip(A) is gauge invariant, and consequently obeys Gauss law. The states 9(A) 
also satisfy Gauss law but they are not gauge invariant. Rather under a gauge 
transformation g, belonging to the homotopy class n g , they change according to 

9(A) -► 9(A 9 ) = e _i9n «*(A) , (88) 

This is the origin of the famous vacuum 6>-angle. 

In 3-dimensional physics, the Chern-Simons term can be added to the ordinary 
Lagrange density of Yang-Mills theories leading to the Lagrange density of the 
topologically massive theories, which is given by 

C m = ^trF^F^ + 8TT 2 kn(A) . (89) 

C m is not gauge invariant; however, at the classical level C m describes a gauge 
invariant theory because the equations of motion that follow from C m are gauge 
invariant: 

V" + ^ a(3 F aP - . (90) 

At the quantum level, the quantity that defines the field theory is the phase ex- 
ponential of the action, and one sees 20 that it is gauge invariant for all values of k 
when the gauge group is Abelian, whereas in the non-Abelian case gauge invariance 
requires that k be quantized as follows: ink = integer. 

Interestingly, one can show that the gauge bosons described by C m have mass k 
(therefore these theories, like the Higgs models, are counterexamples to the naive 
statement that gauge invariance requires gauge bosons to be massless), and this is 
the reason for the name topologically massive theories. 



Q(A) is also related to the Pontryagin density P. In fact, 



P = -^tre^F^Fcf, = d a Cl a (A) , (91) 

where 

n a (A) = -l^e^Hr^AvAp + ^A^Ap) , (92) 

and clearly, for fixed a, fl a (A) can be identified with (±)Q(A) if the dependence of 
the gauge fields on the a-th coordinate is suppressed. 

Numerous mathematically interesting applications have been found for the 
Chern-Simons term; for example, the correlators of Wilson lines in a pure Chern- 
Simons theory (the theory discussed in the following subsection) are related to the 
polynomial invariants of knot theory 21 . 

5.3. Pure Chern-Simons Theory 

The mathematical structure which is important in the context of the hard loops 
high temperature approximation arises in pure Chern-Simons theory. The Lagrange 
density C C s of this theory is obtained as the k — > oo limit of £ m , i.e. by neglecting 
in C m the Yang-Mills kinetic term. Therefore C C s is given by 

C cs = 8ir 2 m(A) . (93) 

The equation of motion which follows from C C s is 

e a » v F„„ = . (94) 

Classically this theory is trivial: the solutions of the equation of motion satisfy 
F^ = and they are therefore pure gauge. However, as a quantum field theory 
this simple model can have some interesting structure. We analyze the theory as a 
canonical quantum field theory in the gauge A a = 0; the Lagrange density is 

Ccs = ^trAiAj . (95) 

By varying with respect to A t one obtains the equations of motion 

Ai = . (96) 

The Hamiltonian vanishes: H = 0. Because of the choice A = 0, the equation gotten 
from Eq.(93) by varying with respect to A a is not obtained; rather a constraint is 
imposed 

G a = ~e ij F% = . (97) 
2 J 

This is the present analog of the Gauss law constraint in ordinary Yang-Mills theory. 
(One easily verifies that the G a 's are the generators of static gauge transformations.) 

There are two ways of quantizing this theory: one can either solve the constraint 
(97) first and then quantize the remaining degrees of freedom or quantize first and 



then solve the constraint. If one solves the constraint first the vector potentials are 
pure gauges, and if the topology of the space is trivial no interesting structure arises 
in the theory. We follow the alternative strategy, i.e. quantize before solving the 
constraint, which does lead to some interesting structures. 

The Lagrange density C C s involves only first order time derivatives; for this 
reason, unlike what happens in Yang-Mills theories, the components of the vector 
potential do not commute: 



[A?(x),A» j (y)] = -e ij 6(x-y) . (98) 

Since H = 0, there is no dynamics, and all the interesting aspects of the theory 
come from the solution of the constraint, i.e. by requiring that the generators G a (x) 
annihilate the physical states 

G a (x)|*>=0. (99) 

We solve Eq.(99) in the Schrodinger picture. Because the spatial components 
of the vector potential do not commute, the wave functionals depend on just one 
of the two spatial components. [Instead for 3-dimensional Yang-Mills theories in 
the (2-dimensional) fixed-time Schrodinger picture, the wave functionals depend on 
both spatial components.] We choose Af, which we call 4> a , as the argument of 
the the wave functionals and A% is realized by functional differentiations with 
respect to 4> a , 

|tf >~*(0) 

A?(x)|* >~0°(x)*(0 (100) 

which is consistent with the commutation relations in Eq.(98). 

In our realization of the Schrodinger picture, the constraint equation (99) can 
be written in the form of functional differential equation 



ft ^ + ^W^- <fcft ^J*^ = ° ■ (101) 

We solve the Gauss law constraint in two steps. First we determine the result 
of a finite gauge transformation g, implemented by the unitary operator U(g), on a 
state, i.e. we evaluate 

U(g)*(4>) = eL Xaix)Ga{x) y{0) = e tG V(4>) , (102) 

[in Eq.(102) we also indicate the explicit form of U(g) in terms of the generators G a , 
and thereby define the operator G] , and then we demand that 

U(g)*{<t>) = e lG *(0) = , (103) 



as required by the Gauss law constraint G a |* >= 0. 



Let us observe that 

G = I X*G« = -j[A'(ft^ + - k I W ° ee G* + 2k I tr&X , (104) 

where is the generator of infinitesimal gauge transformations on the argument 
4> a = A\ of the wave functional. Due to the presence of the term 2k j^tr<j)d 2 \, the 
full action of a gauge transformation on the state is realized with a 1-cocycle (see 
Ref.[22]): 

U{g)^{4>) ee e lG *(0) = e 4G e - lG V G * *(</>) = e lG e~ lG ^((t> 9 ) = e' 2 ^ 011 ^ ^ {(j) 9 ) , (105) 

where (f> 9 = g~ 1 <jyg + g~ 1 dig, and on is obtained by evaluating e iG e - iG <t> with the help of 
the Baker- Hausdorff procedure 

ai (<f>;g) = J tr(2<j>d 2 g g- 1 + g^d^ g- 1 d 2 g) + Ink J u°(g) . (106) 

The uj° which appears in Eq.(106) is the time-like component of the 3- vector uj^, 
which we already encountered in Eq.(86), and one can obtain an explicit expression 
for u>° only after a parametrization for g(x) has been chosen. Alternatively J x ui°(g) 
may be presented as the integral of u(g) = d^(g) over a three-manifold M, whose 
boundary is the two-space on which our Schrodinger functional is defined, 
ai (</>;#) is a one-cocycle because it satisfies the condition 

ai{<t>;g) = ai{(p;gg) - a Y ($ 9 \~g) , (107) 

which is called the one-cocycle condition. 

Eq.(105) concludes the first step in our strategy for solving the Gauss law con- 
straint; at this point we demand that U{g)^{4>) = *(</>), and find that the physical 
wave functionals are not gauge invariant, but rather they satisfy 

= e 2,riai (* ;fl >tf (0) . (108) 

A general wave functional which verifies Eq.(108) can be written in the form 

(0) = e^^VW • (109) 

where ip(4>) is completely gauge invariant [ip{4> 9 ) = ip{4>% W{4>) is given by 

W(<f>) = -k f tr(h- 1 d 1 hh- 1 d 2 h) + 8ir 2 k f w°(h) , (110) 

J x J X 

and h is to be expressed in terms of 4> using 4> = h^dih. W(4>) is a cochain, it satisfies 
W(<t>) = W(4>P)-2ira 1 (<t>;g). 

The term 4irk f x w°(h) is multivalued, and as a consequence it turns out that 
the in order for the wave functional in Eq.(109) to be single valued one needs to 
require that 47rfc=integer. This is the way that the quantization condition, which we 



obtained by requiring the gauge invariance of the phase exponential of the action, 
arises in the present, Schrodinger-picture formalism. 

Note that on a topologically trivial space, one cannot construct a gauge invariant 
functional depending only on 4> a = A\. Consequently ip(4>) in Eq.(109) is a constant 
N, and the general state of Chern-Simons theory is just 

*0) = Ne iWW . (Ill) 

In conclusion, it is interesting to verify that Eq.(97) has indeed been satisfied. To 
this end we consider Ai^ (</>). For i = 1, we have the parametrization (by definition) 
A 1 =<)> = h^dxh. Hence Ai*(0) = ft-^ft*^). For i = 2, A 2 = {l/i)8/8<j) a and A§tf(<£) = 
{8W{4>)/5(j) a )'^{4>). While W is presented in Eq.(llO) as a functional of h, we may 
use the chain rule, and the above formula relating 4> to h. In this way one finds 
A 2 {4>)^{4>) = h^ 1 d 2 h^{4>) and one concludes that, acting on states, A { is a pure gauge 

Ai^((j)) = h^diMlW . (112) 

Therefore, the corresponding field strength vanishes, in keeping with Eq.(97). 

5.4- The Chern-Simons Eikonal 

In 1+1-dimensional point-particle quantum mechanics with Lagrangian 

L = ^f-V(x), (113) 



and energy constraint 



£+n«) = *, aw) 



the eikonal, which is the exponent of the WKB wave functions, is given by 

W(x) = J dx' p(x') , (115) 



where p(x) is the solution of Eq.(114): p(x) = ^2m[E - V(xj\. 
In Chern-Simons theory the Lagrangian can be written as 

L CS = k trA 1 A 2 , (116) 

the Hamiltonian vanishes, and there is the constraint 

F1 2 = d x A a 2 - d 2 Al + f abc A\A c 2 = . (117) 

Clearly there is an analogy with the 1+1-dimensional point-particle quantum me- 
chanics based on the identifications A 1 ~ x and kA 2 ~ p. One can therefore introduce 
WKB wave functional 



m — fif dm* ah*). 



(118) 



where A%(A{) is to be obtained by solving Eq.(117). 

An explicit solution A%(A{) of Eq.(117) is obtained perturbatively in Ref.[19] as 
a series in A\. Here, we observe that W(<f>) coincides with the phase exponential of 
the wave functional that solves the Gauss law constraint (Eqs.(109)-(110)), i.e. the 
WKB/eikonal wave functional is exact in Chern-Simons theory, owing to the simple 
dynamical structure of that theory. To see the result, observe that W(4>) defined by 
Eq.(118) satisfies 

and as a consequence of Eq.(117), W{<j>) also satisfies 

This shows that the wave functional in Eq.(118) satisfies Eq.(lOl). 
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